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REPORT  SUMMARY 


The  research  described  in  this  report  was  carried  out  under  the 
direction  of  Professors  W.  Flygare,  Jr.  and  J.  Jonas  of  the  Department  of 
Chemistry,  Professors  M.  Klein  and  W.  McMillan  of  the  Department  of  Physics, 
and  Professors  N.  Holonyak,  Jr.,  C.  T.  Sah,  and  G.  Stillman  of  the  Depart- 
ment of  Electrical  Engineering  as  an  integral  part  of  the  interdisciplinary 
research  program  of  the  Materials  Research  Laboratory  of  the  University  of 
Illinois . 

Technical  Problem 

The  research  described  in  this  report  investigates  two  areas  of  major 
concern  for  the  development  of  improved  materials. 

(a)  Research  on  Semiconductors:  The  principal  concern  is  to  develop  improved 

light  emitting,  laser,  semiconducting  diodes,  near  infrared  photodiode  de- 
tectors, and  control  defects  and  impurities  in  silicon  diodes,  transitors  and 
integrated  circuits.  Emphasis  has  been  on  the  properties  of  the  compound 
semiconductors  GaAs , InP , GaAs  P,  . In  Ga  P , In,  Ga  P As  for  light 
emitters  and  detectors  and  on  silicon  for  solid  state  circuit  devices. 

(b)  Research  on  Disordered  Materials:  Emphasis  has  been  placed  on  the  de- 

termination of  the  dynamic  structure  of  liquids  including  water,  polymers, 
liquid  crystals,  superionic  conductors,  layer  compounds  and  hydrogen  in  metals 

General  Methodology 

(a)  Research  on  Semiconductors:  Diodes  and  single  and  double  heterojunctions 

were  fabricated  by  liquid  phase  expitaxial  growth.  Optical  absorption  and 
luminescence  measurements  were  made  to  determine  the  performance  of  the 
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devices  and  clarify  the  electronic  mechanisms  involved  in  their  operation. 

A novel  capacitance  technique  was  used  to  determine  the  concentration  and 
properties  of  defects  introduced  into  silicon  devices  during  thermal  pro- 
cessing and  integrated  circuit  manufacture. 

(b)  Research  on  Disordered  Materials:  Sophisticated  light  scattering  tech- 

niques, nuclear  magnetic  resonance,  viscosity  and  elastic  constant  measure- 


ments as  a function  of  pressure  were  the  principal  techniques.  The  experi- 
mental research  was  closely  coupled  with  extensive  theoretical  developments. 


Technical  Results 


(a)  Research  on  Semiconductors 


Laser  diodes  of  InGaPAs  which  operate  in  the  visible  yellow  region  of 
the  spectrum  were  developed.  These  laser  diodes  operate  at  shorter  wave- 
iengths  than  any  other  III-V  compound  system. 

Liquid  phase  expitaxial  layers  of  controlled  thickness  and  properties 
on  GaAs  and  InP  substrates  for  photodiode  detector  devices  were  prepared  and 
the  initial  evaluation  of  the  device  characteristics  begun. 


The  properties  of  tantalum,  gold,  and  aluminum-oxygen  centers  in 
silicon  were  determined. 


(b)  Research  on  Disordered  Materials 

Intermolecular  forces  and  orientational  pair  correlation  of  molecules 
in  liquids  were  determined  and  correlated  with  electric  dipole  moments. 


Molecular  siEe  and  shape  were  found  to  determine  primarily  the  orientational 
pair  correlations. 

The  limitations  of  the  hard  sphere  model  of  a liquid  were  quantitatively 
revealed  by  NMR,  self-diffusion,  viscosity  and  pressure-volume-temperature 


data.  Angular  momentum  correlation  time  data  showed  that  the  effective  hard 
sphere  diameter  of  molecules  is  temperature  dependent.  The  Stokes-Einste in 
equation  was  found  to  be  remarkably  precise  with  the  slipping  boundary  con- 
dition well  satisfied.  The  molecular  conditions  that  yield  these  results 
were  clarified  by  the  data.  At  high  pressure,  the  interesting  result  was 
obtained  that  the  anomalous  properties  of  water  revert  to  more  normal  be- 
havior with  major  changes  in  the  hydrogen  bond  network.  The  effect  of  pres- 
sure on  the  glass  transition  of  several  elastomers  including  natural  rubber 
was  determined. 

Raman  scattering  studies  yielded  detailed  information  concerning  the 
silver  ion  mobility  in  the  superionic  conductor  P-Ag2HgT^  and  AgT.  Informa- 
tion concerning  interatomic  forces  and  bonding  in  titanium  carbide  were  ob- 
tained by  this  technique.  The  growth  of  oxide  layers  on  niobium  was  studied 
by  Raman  scattering  and  scattering  from  hydrogen  in  palladium  observed. 

Proposed  theories  of  the  critical  behavior  near  the  smectic  A to  nematic- 
phase  transition  of  liquid  crystals  were  shown  to  be  faulty.  Theoretical 
molecular  models  of  the  smectic  C and  F phases  were  completed. 

Experimental  data  were  obtained  that  exposed  the  inadequacy  of  current 
theories  of  charge  density  waves  in  the  transition  metal  dichalogenides  and 
A15  compounds.  The  reasons  for  the  failure  of  the  theories  was  determined  to 
be  primarily  the  neglect  of  phonon  entropy.  As  a result,  a satisfactory  theory 
has  been  developed.  A two  dimensional  version  of  Landau's  quantum  hydro- 
dynamics was  used  to  establish  a theory  of  superfluidity  in  two  dimensional 
submonolayer  He‘*  films. 
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Special  Comments 

Examination  of  the  individual  project  reports  which  follow  reveaJ  the 
extensive  pattern  of  collaboration  among  the  investigators  supported  by  this 
contract.  This  collaboration  included  the  preparation  and  sharing  of  mate- 
rials, the  performance  of  measurements,  and  cooperation  in  the  development 
of  theories  based  on  the  experimental  data.  Collaboration  was  not  confined 
to  the  groups  supported  by  this  contract  but  extended  to  important  inter- 
actions with  groups  and  individuals  supported  by  other  Federal  agencies. 

The  semiconductor  research  of  Prof.  N.  Holonyak  benefited  immensely  as  a 
result  of  cooperation  with  investigations  of  the  Monsanto  Chemical  Company 
and  the  liquid  crystals  and  phase  transition  work  similarly  was  assisted  by 
interactions  with  the  Bell  Telephone  Laboratories. 

Of  the  77  publications  resulting  from  the  research  supported  by  the 
contract,  26  publications  or  34%  were  coauthored  by  collaborating  in- 
vestigators . 


Final  Project  Report 

Luminescence,  Lasers,  Carrier  and  Impurity  Effects  in  Compound  Semiconductors 
N.  Holonyak,  Jr. 


Technical  Problem 

In  this  project  we  have  been  concerned  with  two  areas  of  work:  1)  The 

quaternary  In,  Ga  P As  has  been  grown,  via  liquid  phase  epitaxy  (LPE),  and 
fabricated  into  single  and  double  heterojunction  laser  diodes  that  extend 
further  into  the  visible  (yellow)  than  does  any  other  ITI-V  system.  As  part 
of  this  work,  we  have  studied  the  behavior  of  the  quaternary  and  its  two 
ternary  boundaries  GaAs,  P and  In,  Ga  P.  2)  Because  of  its  practical  im- 
portance  (in  making  the  indirect  crystal  quasi-direct)  and  its  fundamentally 
different  behavior  in  a ternary  than  in  a binary,  we  have  studied  further  the 
behavior  of  the  nitrogen  trap  in  GaAs^_^P^;  this  has  been  accomplished  by 
means  of  absorption,  photoluminescence,  and  high  pressure  measurements,  the 
last  on  homo-  and  heteroj unctions . We  have  established  that  the  N trap  in 
GaAs^  yP  , and  in  In^  ^Ga^P,  exhibits  both  short  range  and  long  range  behavior. 


Technical  Results 


InGaPAs  HETEROJUNCTIONS 


Although  it  is  certain  that  In..  Ga  P1  As  can  he  grown  via  vapor 

1“ X X 1*Z  z 

phase  epitaxy  (VPE) , in  this  work  it  has  been  more  convenient  to  use  LPE  to 


grow  single  and  double  heterojunctions  on  GaAs^  ^Py  substrates.  For  both 
single  and  double  heterojunctions  we  have  employed  a cylindrical  slider  boat 


(positioned  in  a small  temperature  gradient)  and  have  determined  that  melt 
wipe-off  after  the  growth  of  each  layer  is  perhaps  the  most  important  LPE 
process  problem,1  particularly  for  In-rich  melts  suitable  for  large  energy 


gap  layers. 
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For  the  case  of  single  heterojunctions,  the  active  region  of  the 
structure  lies  in  the  GaAs,  P substrate . 2 * 3 This  has  made  it  possible 

l-y  y 

to  study  a number  of  the  fundamental  properties  of  the  N trap  in 

GaAs,  P in  the  region  near  the  direct-indirect  transition. 4-8  We  mention 
l-y  y 

that  as  part  of  this  work  we  have  determined,  via  the  first  successful 

direct-gap  exciton  absorption  measurements  in  a I1I-V  alloy,  *»10  that  the 

direct-indirect  transition  in  GaAs,  P lies  at  v = y„  0.45  (77°K)  and 

l-y  y c 

in  In,  Ga  P at  x = x ^ 0.72  (77°K) . 

1-x  x c 

Utilizing  GaAs^  ^P^  as  a substrate,  we  have  succeeded  in  growing 

In,  Ga  P As  DH  laser  diodes  that  operate  in  the  yellow  at  77°K.11>1? 

1 “X  X l— z z 

These  laser  diodes  have  been  operated  also  in  external  grating  cavities 
and  have  demonstrated  to  what  extent  ''homogeneous11  or  "inhomogeneous"  broad- 
ening can  be  expected  in  a diode  laser.13*14  As  the  LPE  crystal  growth  art 
has  improved,  it  has  proven  possible  to  operate  In^^Ga^P^^As^  OH  laser 
diodes  (vl40  meV  heterobarrier)  at  room  temperature  with  pulsed  current 
densities  of  '^2xl04A/cm2 . 1 5 The  heterobarrier  height  can  be  increased  to 
250  meV  without  sacrificing  the  capability  of  laser  operation  (300°K)  at 

O i r 

6500  A.  c These  results  indicate  that  when  more  sophisticated  stripe- 

geometry  diode  structures  are  employed,  and  the  LPE  layer  growth  is  further 

developed,  it  should  be  possible  to  operate  red-orange  In  Ga  P As  DH 

1“X  X I— z z 

laser  diodes  continuously  at  room  temperature.  The  composition  (xc,  z ) at 
which  these  quaternary  laser  diodes  become  limited  by  the  direct-indirect 
transition  (xc-0.52zc  = 0.73,  77°K;  xc-0.51zc  = 0.72,  300°K)  has  been  de- 
termined by  pressure  measurements17  and  by  observation  of  diode  "freeze-out" 

(300+77°K)  or  impedance  increase  for  lattice-matched  wide-gap  In,  Ga  P,  As 

1-x  x 1-z  z 

(x  "•  xc,  z f-,  0.01)  layers  grown  on  successively  higher  composition  GaAs^_^P^ 
substrates . 1 


w 


I 
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NITROGEN  TRAP  IN  III-V  ALLOYS 

By  means  of  absorption,  photoluminescence  and  electroluminescence 

measurements  on  In,  Ga  P:N10>ig>20  and  by  a large  series  of  measurements 
1— x x 

on  GaAs^_^Py : N,  7»  R»  2 1-/ 6 we  have  shown  that  the  previous  assignment  of  an 
NN-pair  recombination  transition27  in  these  alloys  (analogous  to  GaP:N  be- 
havior) is  in  error,  and  that  this  transition,  now  labeled  N , is  due  to 

A 

recombination  at  a single  N site.  Theory,  with  measurement,8  shows  that  the 

N trap  in  GaAs,  P (and  in  In,  Ga  P)  exhibits  a combination  of  long-range 
1-y  y 1— x x 

and  short-range  behavior,  the  latter  being  most  important  for  light  emitting 

purposes  (indirect  crystal).26’27  The  two  long-range  states,  N and  N , are 

I A 

associated  with  the  T and  X band  minima.7’8  The  state  N , which  is  useful 

in  indirect  crystal  for  light  emitting  diodes  (LED's),2  is  basically  fairly 

weak  except  for  the  oscillator  strength  coupled  into  this  transition  by  the 

short-range  state,  N.fi>28  Pressure  measurements  on  In,  Ga  P,  As 

1-x  x 1-z  z 

GaAs^  yPy'N  heteroj unctions2 » 7 > 2 8 in  the  composition  range  where  the  short- 

range  trap,  N,  is  degenerate  with  N and  where  N is  degenerate  with  N in- 

I A I 

dicate  strong  interaction  and  splitting  of  these  states  at  the  crossover 
compositions . 7 


Implications  for  Future  Research 

We  have  demonstrated  that  the  quaternary  In,  Ga  P,  As  can  be  used 

J 1-x  X 1-z  z 

to  fabricate  single  and  double  heterojunctions,  the  latter  able  to  operate 
well  into  the  yellow  portion  of  the  spectrum.  This  is  the  highest  energy 
diode  laser  system  yet  demonstrated.  Although  it  is  not  yet  proven,  the 
In.  Ga  P.  As  system  potentially  is  capable  of  lattice-latching  when  one 
layer  is  grown  epitaxially  on  another.  For  technological  reasons  it  is  im- 


portant that  this  alloy  is  free  of  A1  and  its  difficult  oxides. 
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Based  on  experiments  conducted  on  In,  Ga  P and  GaAs,  P , we  have 
v 1-x  x 1-y  y 

been  able  to  show  that  the  N trap  in  these  alloys  exhibits  a combination 
of  long  range  and  short  range  behavior  and  that  NN-pair  recombination 
transitions  (unlike  in  GaP:N)  are  not  important.  This  work  is  useful  in 
helping  to  clarify  the  recombination  mechanisms  that  are  important  in 
certain  classes  of  LED's. 

Both  of  these  areas  of  study  are  far  from  complete  and  require  more 
work.  Based  on  the  work  reported  here,  it  is  likely  that  a CW,  room  tern- 

O 

perature  visible  (\-6500A)  DH  diode  laser  can  be  built;  also,  it  is  likely 
that  the  work  we  describe  here,  pursued  further,  will  improve  LED's  by  a 
factor  of  10. 
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Final  Project  Report 

Compound  Semiconductors  for  Near  Infrared  Detectors 
G.  Stillman 

Technical  Problem 

This  project  was  initiated  in  July  1975,  to  study  various  III-V  com- 
pounds for  applications  in  near  infrared  photodiode  detectors.  The  performance 
of  silicon  photodetectors  that  are  available  commercially  is  limited  bv  a 
response  time-quantum  efficiency  product  at  wavelengths  longer  than  about 
0.9  pm  because  of  the  low  absorption  coefficient  of  silicon  in  this  wave- 
length range.  Photodiodes  in  direct  band  gap  III-V  compound  semiconductors 
do  not  have  the  same  limitation  and  are  inherently  capable  of  much  shorter 
(transit  time  limited)  response  times  at  a given  quantum  efficiency.  For 
highest  sensitivity  in  wide  bandwidth  detection  applications  it  is  also  de- 
sirable to  operate  the  detectors  with  avalanche  gain,  and  the  larger  the  dif- 
ference in  electron  and  hole  ionization  coefficients,  the  higher  the  sensi- 
tivity that  can  be  obtained  in  a properly  designed  device.  The  initial 
objective  of  this  project  was  to  establish  a laboratory  for  the  liquid  phase 
epitaxial  growth  of  GaAs  and  InP,  the  general  characterization  and  evaluation 
of  semiconductor  material,  and  for  the  fabrication  and  evaluation  of  photo- 
diode detectors.  The  second  objective  was  to  fabricate  special  device  struc- 
tures which  will  permit  the  accurate  determination  of  the  electron  and  hole 
impact  ionization  coefficients  over  a wide  range  of  electric  fields  above  and 
below  room  temperature.  Besides  being  of  interest  for  the  design  or  avalanche 
photodiodes,  these  parameters  are  also  of  great  importance  in  IMPATT  device 
analyses . 


General  Methodology 


The  slider  boat  liquid  phase  epitaxial  technique  was  chosen  for  the 
initial  growth  of  GaAs  and  InP  to  be  used  in  near  infrared  detector  applica- 
tions because  of  its  versatility  for  multiple  layer  growth  and  because  of 
its  relative  simplicity.  Separate  systems  have  been  set  up  for  the  growtli 
of  GaAs  and  InP.  The  InP  LPE  svstem  incorporates  a transparent  gold  furnace 
to  permit  accurate  determination  of  the  liquidus  and  solidus  temperatures 
for  a given  melt  composition.  The  growth  techniques  presentlv  being  utilized 
include  supersaturation,  constant-temperature,  and  ramp  down  growth. 

For  the  characterization  of  the  semiconductor  material  grown,  facil- 
ities have  been  set  up  for  Hall  effect  measurements  at  room  temperature  and 
liquid  nitrogen  temperature,  and  for  C-V  measurements  and  analyses  using 
Schottkv  barrier  diodes.  In  addition,  optical  facilities  exist  for  the 
measurement  of  transmission  and/or  absorption  over  the  wavelength  range  from 
the  near  ultraviolet  out  to  about  50  pm,  and  these  measurement  techniques 
are  used  to  evaluate  the  quality  of  the  material  (deep  levels,  etc.)  and  to 
determine  the  band  gap  of  semiconductor  alloys. 

The  fabrication  of  wide  bandwidth  photodiode  detectors  requires  the 
use  of  photolithographic  processing  techniques,  and  a small  facility  which 
permits  this  processing  has  been  set  up  in  collaboration  with  Professor 
Nick  Holonyak,  Jr.  Included  in  this  facility  are  photoresist  and  two  inch 
mask  alignment  equipment,  and  the  capability  of  depositing  pyrolytic  SiO^, 
SiON  and  Si^N^  for  use  in  passivation,  pattern  definition,  and  masking  for 
diffusion  and  etching.  An  alloy  station  which  uses  an  HCI-H2  atmosphere  is 
employed  for  the  low  temperature  preparation  of  ohmic  semiconductor  contacts. 

The  characterization  of  the  photodiode  detectors  includes  measurements 
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of  quantum  efficiency,  relative  spectral  response,  noise,  and  response  time. 
For  the  measurement  of  quantum  efficiency,  a spot  measurement  is  made  using 
a HeNe  laser  and  a calibrated  photodiode.  The  spectral  quantum  efficiency 
and  the  spectral  responsivity  are  then  determined  by  using  the  spot  measure- 
ment to  calibrate  relative  spectral  measurements  made  using  a thermocouple 
reference  detector  and  a monochrometer  with  a tungsten  light  source.  Standard 
laboratory  techniques  are  used  for  the  noise  and  response  time  measurements. 

Technical  Results 

Liquid  phase  epitaxial  layers  with  good  morphological  and  electrical 
characteristics  can  now  be  grown  routinely  on  <100>  and  < 1 1 1 > N+  and  high 
resistivity  Cr-doped  substrates.  There  is  generally  good  agreement  between 
the  experimental  layer  thickness  and  that  calculated  using  diffusion  limited 
growth  models,  so  that  the  layer  thickness  can  be  controlled  quite  accurately 
for  thicknesses  greater  than  1 um.  This  is  important  for  the  design  and 
fabrication  of  both  photodiodes  and  particular  structures  for  electron  and 
hole  ionization  coefficient  measurements.  The  doping  can  generally  be  con- 
trolled quite  well  in  the  low  1015  cm-3  range,  as  confirmed  by  Hall  coeffi- 
cient measurements  of  the  layers  grown  on  insulating  substrates  and  by  C-V 
measurements  of  the  layers  on  conducting  substrates.  This  growth  technique 
is  now  being  used  to  fabricate  special  structures  for  the  measurement  of 
electron  and  hole  impact  ionization  coefficients.  As  the  InP  growth  param- 
eters are  determined,  these  measurements  will  be  extended  to  this  material , 
and  following  this  to  InGaAsP  quaternary  epitaxial  layers  for  longer  wave- 
length near  infrared  applications. 
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Final  Project  Report 

Properties  of  Recombination  Centers  in  Semiconductors 

a 

C.  T.  Sah 

, 

Technical  Problem 

This  project  studies  the  recombination  and  diffusion  kinetic  properties 
of  impurity  and  defect  recombination  centers  in  silicon  introduced  during 
crystal  growth  and  subsequent  high  temperature  processing  steps  commonly  used 
in  integrated  circuit  manufacturing.  The  capacitance-voltage  and  transient 
capacitance  techniques  are  used  to  determine  thermal  emission  rates  of  elec- 
trons and  holes,  the  energy  level  and  the  concentration  profile  of  these 
centers.  Detailed  studies  were  made  for  Ta  and  aluminum-oxygen  centers  as 
well  as  low  temperature  diffusion  of  gold  into  an  aluminum  sink  in  silicon. 


Technical  Results 

Detailed  measurements  of  the  parameters  characterizing  the  generation 
and  recombination  properties  of  electrons  and  holes  at  Ta  centers  in  silicon 
have  been  made  and  reported.1  High  concentrations  of  Ta  were  introduced  into 
the  silicon  melt  but  only  a small  fraction  is  electrically  active  in  the 
grown  silicon  crystal. 

The  effect  of  the  stationary  trapped  charges  at  recombination  centers 
located  at  the  edge  of  the  space  charge  region  of  p-n  or  Schottky  barrier 
junction  has  been  analyzed  and  determined  experimentally.  The  model  including 


this  stationary  charge  has  been  used  to  evaluate  the  spatial  variation  of 
concentration  of  the  recombination  centers  in  the  region  very  near  the  p-n 
junction  or  the  metal-semiconductor  interface.  One  interesting  result  ob- 


the 


tained  is  the  high  diffusion  rate  of  gold  from  silicon  into  the  aluminum  sink 
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at  450°C  in  an  A1  on  n-Si  Schottky  barrier  diode  structure.  It  shows  that 
gold  outdif fuses  into  the  aluminum  sink  following  the  simple  diffusion  law 
which  gives  an  erf  concentration  profile  with  dif f usivities  consistent  with 
those  extrapolated  from  800  to  1200°C.2 

The  capacitance-voltage  method  has  been  applied  to  the  determination 
of  impurity  concentration,  resistivity  and  conductivity  mobility  in  p-type 
silicon.  The  impurity  and  carrier  concentrations  are  obtained  from  the 
capacitance-voltage  measurements  made  on  Al  on  n-Si  Schottky  barrier  diodes 
while  the  resistivity  is  measured  on  an  ohmic  diode  adjacent  to  the  Schottky 
diode.  The  analyses3  of  the  data  had  neglected  the  deionization  effect  al 
low  temperatures  which  is  not  justifiable.  A revised  calculation  of  the 
experiment  has  been  undertaken  and  will  be  submitted  for  publication. 

The  capacitance  transient  technique  has  also  been  applied  to  the  study 
of  the  traps  produced  during  heat  treatment  in  aluminum-doped  silicon.  There 
has  been  some  interest  in  Al-doped  p-type  silicon  in  certain  IR  applications 
due  to  the  larger  thermal  activation  energy  of  aluminum  in  silicon.  Detailed 
measurements  have  shown  that  the  defect  centers  detected  are  aluminum-oxygen 
complexes  which  produce  three  hole-trap  levels  and  one  electron-trap  level  in 
the  silicon  band  gap.  The  complexes  break  up  at  1000°C  or  above  and  the 
oxygen  member  of  the  complex  diffuses  toward  a heavily  phosphorus  doped  layer. 
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Final  Project  Report 

Liquid  Crystals  and  Phase  Transitions 
W.  McMillan 


Technical  Problem 

This  research  was  concentrated  in  two  areas:  liquid  crystals  and 

charge  density  waves. 

In  liquid  crystals  the  technical  problem  is  to  gain  a microscopic 
understanding  of  the  liquid  crystal  phases  and  their  phase  transitions. 

The  principal  experimental  method  used  is  Rayleigh  light  scattering  and 
a new  optical  method  for  measuring  elastic  constants.  The  principal 
theoretical  tools  have  been  microscopic  mean  field  theory,  phenomenologi- 
cal Landau  theory,  and  renormalization  group  methods.  We  have  perfected 
a new  high  precision  optical  method  for  measuring  elastic  constants  and 
viscosity  coefficients  of  liquid  crystals  in  the  nematic  phase.  We  have 

used  this  method  to  study  the  critical  behaviour  near  the  smectic  A to 

nematic  phase  transition  in  CBOOA  (cyanobenzy lidene-octv loxyani 1 ine) . 

These  new  results,  together  with  our  Rayleigh  scattering  measurements, 
show  that  all  the  proposed  theories  of  the  critical  behaviour  are  incor- 
rect. On  the  theoretical  side  we  have  completed  a series  of  papers  on  the 

Landau  theory  approach  on  several  phase  transition  problems. 

General  Methodology 

In  the  study  of  charge  density  waves  (CDW)  in  metals  the  thrust  has 
been  mainly  theoretical,  in  collaboration  with  experimenters  Fred  Brown, 
Bob  Graven,  Miles  Klein  and  the  Bell  Labs  group.  The  overall  objective  is 
to  construct  a quantitative  theory  of  CDWs  in  the  transition  metal 
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dichalcogenides  (2H-TaSe^ ) and  to  test  it  by  comparison  with  experiment. 

A Landau  approach  has  been  very  useful  in  describing  the  phase  transitions 
and  the  defects  (CDW  dislocations  and  discommensurations) . We  have  re- 
cently shown  that  previous  microscopic  theories  cannot  describe  CDWs  in 
2H-TaSe9  and  have  reformulated  the  theory.  We  have  calculated  lattice 
dynamics  in  the  distorted  phase  in  order  to  interpret  Klein's  Raman  data. 

In  addition  we  have  studied  superlattice  formation  in  lT-TiSe2  with  x rays. 

Technical  Results 

We  measured  the  twist  elastic  constant  and  the  viscosity  of 

CBOOA  near  the  smectic  A-nematic  transition  using  the  standard  Rayleigh 
scattering  technique  and  with  temperature  control  and  homogeneity  better 
than  .001°C  (at  83°C).  Our  results  fit 


K22(T)  - K2;  + C/(T  - Tan) 


with  v = 0.5  with  a similar  expression  for  y . This  result  contradicts 
deGennes  analogy  to  liquid  helium  (v  = 2/3)  and  Halperin  et  al.'s  predic- 
tion of  a first  order  transition. 

We  have  developed  a new  experimental  technique  of  forcing  a small 
distortion  on  the  liquid  crystal  using  an  external  magnetic  field  and 
measuring  the  distortion  using  a sensitive  optical  interference  technique. 
The  method  yields  higher  precision  than  Rayleigh  scattering  and  is  useful 
within  .003°C  of  the  phase  transition.  The  standard  Freedericksz  method 
imposes  large  distortions  and  fails  .1°C  above  the  transition.  Our  results 
for  the  band  elastic  constant  are  similar  to  (1)  with  v = .73.  This  re- 
sult contradicts  mean  field  theory  (v  = .5).  We  have  therefore  shown  that 
all  the  theories  of  the  critical  behaviour  are  wrong.  These  results  also 
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indicate  that  there  is  an  essential  anisotropy  in  the  critical  behaviour 
with  the  transverse  and  longitudinal  correlation  length  having  different 
exponents,  a result  suggested  by  ray  early  x-ray  study  and  by  Litster's 
light  scattering  study  in  the  smectic  A phase.  The  experiments  indicate 
that  the  theory  will  have  to  be  reformulated  keeping  the  anisotropy  in  mind. 

The  molecular  theory  of  the  smectic  phases  is  essentially  complete; 
a model  of  the  smectic  F phase  will  soon  be  published.  The  molecular  model 
of  the  smectic  C phase  (orientational  ordering  due  to  dipole-dipole  inter- 
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actions)  is  still  controversial.  However,  experiments  in  progress  (Doane 
at  Kent  State  and  Blinc  in  Yugoslavia)  should  settle  the  question. 

In  a series  of  papers  we  have  established  the  Landau  theory  approach 
to  charge  density  waves  in  transition  metal  dichalcogenides  and  in  A 15 
compounds.  This  approach  has  been  very  useful  and  very  successful.  With 
the  completion  of  the  heat  capacity  (Cravens)  and  Raman  (Klein)  experiments 
in  Urbana  there  is  enough  experimental  data  to  test  the  proposed  micro- 
scopic theories  of  CDWs.  The  proposed  theories  all  fail  to  describe  the 
data.  The  technical  reason  that  the  theories  fail  is  that  they  assume  a 
long  correlation  length  in  a BCS-like  energy  gap  model.  The  theories  simply 
neglect  the  phonon  entropy.  We  have  shown  that  for  2H-TaSe2  the  phonon  en- 
tropy is  dominant  and  have  reformulated  the  theory  accordingly.  The  next 
step  is  to  include  the  correct  band  structure  in  the  theory  and  to  check  it 
against  the  dellaas-van  Alphen  (Graebner)  and  deHaas-Shubnikov  (Coleman) 
data.  At  least  for  the  dichalcogenides  the  experiments  seem  to  be  pushing 
us  rapidly  toward  a viable  microscopic  theory.  A very  fruitful  collabora- 
tion with  Miles  Klein  on  the  Raman  experiments,  Bob  Craven  on  the  heat  ca- 
pacity measurements  and  Fred  Brown  on  the  study  of  lT-TiSe^  has  developed 
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within  the  Materials  Research  Laboratory. 

We  have  made  a great  deal  of  progress  in  understanding  superfluidity 
in  submonolaver  He4  films  in  the  last  year  due  to  a collaboration  with 
Jack  Mochel.  Mochel  has  very  beautiful,  high  precision  measurements  of 
third  sound  velocity  and  attenuation  from  . 1°K.  to  1°K.  We  have  interpreted 
these  measurements  using  a theoretical  model  whicli  is  a two-dimensional 
version  of  Landau's  quantum  hydrodynamics.  The  model  predicts  an  elementary 
excitation  spectrum  which  can  be  used  to  calculate  the  temperature  dependence 
of  third  sound,  thermal  boundary  resistance  and  critical  velocity  effects. 
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Final  Project  Report 

Light  Scattering  from  Non-Vlbrat ional  Modes  in  Solids 
M.  S.  Klein 

Technical  Problem 

The  main  objective  of  this  work  has  been  to  study  non-vibrational 
motions  in  disordered  solids.  Particular  emphasis  has  been  placed  on  two 
classes  of  system  in  which  atoms  or  ions  diffuse  very  rapidly:  hydrogen 

in  metals  and  fast-ion  conductors.  The  technique  used  is  Raman  scattering, 
which  has  not  been  applied  extensively  to  metals.  A secondary  objective  of 
this  work  has  been  to  broaden  knowledge  of  Raman  scattering  in  metals  by 
applying  the  technique  to  study  lattice  vibrations  in  disordered  metals  and 
in  metals  that  have  electronically  driven  phase  transitions. 

C’.eneral  Methodology 

Our  general  technique  is  laser-excited  Raman  scattering.  Because  of 
the  great  deal  of  disorder  in  the  materials  under  investigation,  much  of  the 
work  requires  a third  monochromator  to  be  added  to  the  usual  double  mono- 
chromator to  cause  a further  reduction  in  the  parasitic  elastically  scattered 
laser  light.  Much  of  our  effort  on  hydrogen  in  metaTs  has  involved  the  bcc 
metals  niobium  and  tantalum,  which  oxidize  readily,  especially  at  the  tem- 
peratures necessary  to  obtain  high  equilibrium  hydrogen  concentrations.  Thus 
the  experiments  have  been  carried  out  in  an  ultra-high  vacuum  system.  This 
has  allowed  the  beginnings  of  a study  of  Raman  scattering  from  a thin  oxide 
laver  on  niobium  and  from  adsorbed  molecules  on  the  oxide. 

None  of  the  samples  studied  in  this  work  has  been  readily  available, 
and  a great  deal  of  effort  has  been  expended  in  sample  preparation.  The  fast 
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ion  conductors  studied,  RbAg^I,.  and  NH^Ag^I^,  have  been  grown  from  aqueous 
solution  in  a joint  effort  involving  members  of  three  research  groups: 

Salamon,  Klein,  and  Lazarus.  The  niobium: hydrogen  and  tantalum:hydrogen 
work  has  been  done  in  collaboration  with  Prof.  H.  Birnbaum,  who  has  pre- 
pared the  samples  and  has  been  closely  involved  in  various  treatments  given 
the  samples  during  the  Raman  measurements.  The  crystals  for  the  TiC  measure- 
ments were  supplied  by  Prof.  W.  Williams,  who  has  highly  characterized  them. 

The  samples  of  2H-TaSe2  came  from  Dr.  S.  F.  Meyer  of  Prof.  C.  Slichter's 
group.  Transport  and  specific  heat  measurements  on  them  were  made  by  Dr.  R. 
Craven  of  IBM,  formerly  a member  of  Dr.  Salamon 's  group.  The  samples  of 
lT-TiSe0  were  prepared  by  K.  C.  Woo  of  Dr.  F.  Brown's  group.  The  Raman  and 
infrared  work  on  TCNQ  was  performed  on  samples  prepared  by  Prof.  G.  DePasquali. 

Technical  Results 

(a)  Pre-resonance  Raman  spectra  have  been  obtained  by  Professor  Klein 
working  with  Dr.  Wozniak  and  Professor  DePasquali  for  TCNQ  and  LiTCNQ  in 
acetonitrile  solution  using  an  Ar+  - Kr+  laser  and  a tunable  rhodamine  6G 
dve  laser.'  The  theory  of  Albrecht  and  Hutley  has  been  used  to  calculate  fre- 
quency factors  for  the  intensity  variations  for  several  symmetric  vibrational 
modes  of  each  molecule.  The  observed  spectra  for  TCNQ  and  LiTCNQ  with  violet, 
blue,  and  green  excitation  give  evidence  of  resonance  enhancement  due  to 
vibronic  mixing  between  at  least  two  violet  and  ultraviolet  transitions.  The 
Raman  spectra  for  LiTCNQ  with  yellow,  orange,  and  red  excitation  show  enhance- 
ment due  to  a single  electronic  excitation  in  the  near  infrared. 

(b)  Professor  M.  Klein  has  developed  a theory  of  light  scattering  from 
the  relaxation  modes  of  rapidly  diffusing  atoms  or  ions  using  the  independent 
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Professor  Klein  and  D.  Gallagher  have  obtained  Raman  data  on  the 
superionic  conductor  Ag^Rbl,.  at  temperatures  of  136  K to  448  K. ‘ At  209  K 
there  is  a second  order  phase  transition  from  a highly  ionic  conducting 
phase  to  one  with  less  conductivity.  The  Raman  spectrum  shows  no  abrupt 
changes  at  the  transition.  At  high  temperatures  there  are  two  main  features — 
a symmetric  peak  centered  at  105  cm-1,  and  a shoulder  on  the  Rayleigh  line 
at  about  23  cm-'.  As  the  temperature  is  lowered,  the  shoulder  resolves  it- 
self into  three  peaks.  Polarized  Raman  spectra  cannot  he  taken  below  the 
transition  at  209  K because  the  crystal  becomes  a random  superposition  of 
uniaxial  ferroelastic  domains. 

The  Ag  ions  are  surrounded  by  a distorted  tetrahedral  cage  of  I ions. 

A similar  local  geometry  is  found  in  8 - Ag2HgI^4  and  in  phase  IV  of  Agl. 

The  Raman  spectrum  of  these  materials  shows  a peak  at  20  - 30  cm-1  and  one 
at  about  110  cm-1.  Tentative  explanations  of  the  spectra  in  all  materials 
are  the  same — the  100  cm-1  peak  is  a "breathing"  node  of  the  I-tetrahedra , 
and  the  20  - 30  cm-1  mode  is  the  "attempt  frequency"  of  the  silver  ion  os- 
cillating in  its  potential  well.  This  is  presumably  the  same  attempt  vibra- 
tion involved  in  the  diffusive  jumping  of  the  Ag  ions  from  site  to  site.  The 
100  cm-’  breathing  mode  is  presumably  the  polaron  mode  responsible  for  dress- 
ing the  Ag  ion's  motion.  Thus  the  two  Raman  features  may  directly  relate  to 
two  processes  that  are  very  important  to  an  understanding  of  the  microscopic 
dynamics  of  this  superionic  conductor. 

(c)  Raman  spectra  have  been  obtained  by  Professor  Klein,  Mr.  Holy,  and 
Professor  Williams  on  samples  of  TiC  containing  various  concentrations  of 
carbon  vacancies.  The  vacancies  "induce"  Raman  scattering  that  would  other- 
wise not  be  allowed.  The  first-order  impurity-induced  spectra  show  broad 
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peaks  at  285  and  410  cm-1  in  the  acoustical  phonon  region,  a gap,  and  peaks 
at  570  and  680  cm  1 . Some  of  the  details  of  the  spectra  depend  on  the 
vacancy  concentration;  the  general  features  do  not.  These  data  reveal  new 
information  about  interatomic  forces  and  bonding  in  these  materials.  For 
instance,  the  data  suggest  strongly  that  the  first  shell  of  C atoms  around 
a C vacancy  (next-nearest  neighbors)  have  their  electronic  polarizability 
more  strongly  affected  by  the  vacancy  than  do  the  first  nearest  neighbor 
Ti  atoms. 

(d)  Raman  spectra  have  been  obtained  by  Professor  Klein  and  Mr.  Holy 
on  IT  and  2H  polytypes  of  TaS2  obtained  from  Dr.  S.  Meyer  of  Professor 
Slichter's  group.  This  material  is  a nearly  2-dimensional  layered  compound 
with  very  weak  bonding  between  layers.  In  the  IT  polytype,  two  Raman  modes 
have  been  seen,  as  predicted  by  symmetry  arguments.  Two  weaker  "forbidden" 
modes  are  also  observed.  These  may  be  caused  by  a distortion  of  symmetry 
breaking  associated  with  a charge  density  wave  known  to  be  present  in 
lT-TaS2  at  room  temperature.  In  2H-TaS2  four  Raman  are  observed;  the  lowest 
one  at  24  cm  1 is  a direct  measure  of  the  interlayer  force  constants. 

Polarized  Raman  spectra  have  been  obtained  on  2H-TaSe2  at  temperatures 
from  300  K to  27  K,  below  the  transition  to  a commensurate  charge-density 
wave  and  superlattice. 6 At  room  temperature  all  the  allowed  Raman  modes  of 
the  2H  structure  are  seen  along  with  a very  strong,  broad  band  at  148  cm- 1 
that  is  assigned  to  a combination  of  two  LA  phonons  having  wave  vectors  on 
the  line  I from  about  0.40  to  1.00  q^.  The  dispersion  curve  for  this  branch 
shows  softening  and  a Kohn  anomaly.7 

Diffraction  studies  on  2H-TaSe27,fi  show  a transition  from  a normal 


lattice  to  one  with  an  incommensurate  charge  density  wave  (CDW)  at  122  K 


followed  at  90  K by  a transition  to  a commensurate  state.  These  studies 
strongly  suggest,  but  do  not  prove,  the  coherent  superposition  of  three 
CDW's.  Neutron  diffraction  studies  give  the  static  displacements  in  the 
distorted  state,  but  they  are  relatively  insensitive  to  the  absolute  phase 
angle.  At  low  temperatures  six  new  low  frequency  Raman  modes  are  seen. 

These  are  interpreted  as  new  zone  center  optical  phonons  that  result  from 
folding  in  the  anomalous  LA  mode  at  2/3q^.  Their  number  and  Raman  sym- 
metries are  consistent  only  with  an  inversion-symmetric,  3CDW  state.  The 
Raman  data  have  removed  the  ambiguities  remaining  from  the  diffraction 
studies  and  allowed  a determination  of  the  space  group  of  the  commensurate 
CDW  state.  W.  L.  McMillan  has  produced  a microscopic  model  of  the  lattice 
dynamics  of  this  state.'  The  parameters  in  this  model  are  fitted  to  the 
frequencies  of  the  four  strongest  new  Raman  lines,  which  are  assigned  to 
amplitude  and  phase  modes  of  the  charge  density  wave. 

Raman  studies  have  been  performed  on  lT-TiSe2  by  J.  A.  Holy  and  Prof. 
M.  Klein  in  collaboration  with  K.  C.  Woo  and  Prof.  F.  Brown.  Below  the 
transition  to  a CDW  state  new  Raman-active  and  infrared-active  phonon  modes 
are  seen.  This  is  a more  complicated  system  than  ZH-TaSe^  because  the  new 
modes  are  derived  from  the  folding  of  more  than  one  phonon  branch.  Detailed 
understanding  may  require  more  far  infrared  measurements  and  the  results  of 
inelastic  neutron  measurements. 

(e)  Attempts  extending  over  a two-year  period  to  see  Raman  scattering 
from  the  vibrations  of  hydrogen  and  deuterium  in  niobium  and  tantalum  have 
not  been  successful.  Raman  scattering  was  seen  from  pure  Ta  in  the  form  of 
nice  two-phonon  density  of  states.  This  altered  somewhat  upon  hvdrogenat ion 
but  no  high  frequency  hydrogen  modes  appeared.  In  the  case  of  tantalum  this 
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negative  result  has  a positive  significance:  the  vibrations  of  the  hydro- 

gen atoms  do  not  appreciably  couple  to  the  optical  susceptibility  at  the 
laser  photon  energy,  2.4  eV.  This  implies  that  the  hydrogen  does  not  in- 
troduce an  electronic  level  within  2.5  eV  of  the  Fermi  level  that  can  be 
reached  from  the  Fermi  level  via  electric  dipole  transitions.  There  was 
some  oxide  contamination  during  the  experiment  on  hydrogen  in  niobium  and 
no  two-phonon  spectrum  was  seen,  thus  the  negative  result  on  that  system  is 
not  as  significant  as  it  is  for  tantalum. 

The  Raman  work  on  niobium  showed  one  broad  peak  at  640  cm  that  is 
probably  due  to  an  oxide  layer.  To  better  understand  this  peak  the  oxida- 
tion of  clean  Nb  surfaces  at  very  low  partial  pressure  of  0^  was  studied. 

At  200°C  an  amorphous  NbO  layer  initially  forms  giving  a very  broad  Raman 
peak.  As  the  oxide  film  thickens  the  NbO  crystallizes  and  the  broad  peak 
develops  into  several  sharp  Raman  peaks  characteristic  of  NbO.  Continued 
exposure  results  in  the  growth  of  an  amorphous  layer  of  Nl^O,.  on  toP  °f  the 
NbO,  and  this  results  in  a broad  Raman  peak.  The  Nb20,_  develops  into  a 
crvstalline  form  as  it  thickens,  and  sharp  Raman  lines  characteristic  of 
Nb2<>^  are  then  observed. 

Higher  frequency  Raman  peaks  were  also  observed  at  1350  and  1600  cm-1 
on  Nb  samples  held  in  a high  vacuum  system.  They  were  not  due  to  intentional 
additions  of  hydrogen.  Their  positions  coincide  exactly  with  some  infrared 
lines  observed  with  deposited  Al  and  Cu  films  exposed  to  formic  acid  vapors.'1 
The  adsorbed  ions  were  thought  to  be  HCC^-.  The  niobium  samples  were  electro- 
polished  in  lactic  acid  before  being  placed  in  the  vacuum  system,  and  lactic 
acid  contains  HCC^  units.  If  such  units  are  responsible  for  the  Raman  lines, 


they  are  probably  formed  on  an  oxide  layer  rather  than  on  the  niobium  itself. 


34 


Implications  for  Further  Research 

There  are  several  immediate  implications  of  the  above  results  that  are 
affecting  the  continuation  of  this  work.  The  Raman  work  on  the  fast  ion  con- 
ductor RbAg^Ij  is  currently  being  extended  to  NH^Ag^I^,  which  has  the  same 
properties  of  ionic  conduction  and  low  frequency  vibrational  modes  but  has  in 
addition  the  internal  vibrational  modes  of  the  ammonium  ion  that  will  be 
affected  by  the  static  and  dynamical  disorder. 

Theoretical  work  has  been  started  and  additional  experiments  are  planned 
to  explain  why  the  charge-density-waves  in  compounds  such  as  TaSe2  produce  such 
strong  new  Raman-active  modes  in  the  low  temperature  phases. 

Attempts  to  observe  Raman  scattering  from  tiydrogen  in  metals  are  being 
extended  to  fee  metals.  A positive  result  has  recently  been  observed  in 
Pd:H.  Additional  feasability  studies  are  planned  to  determine  if  Raman  scat- 
tering can  be  used  to  study  chemisorbed  species  and  thin  oxide  layers  on 
metal  surfaces  under  controlled  conditions. 
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Final  Project  Report 

Higli  Pressure  Nuclear  Magnetic  Resonance  and  Raman  Study  of 
the  Dynamic  Structure  of  Liquids,  Disordered  Solids  and  Polymers 

J.  Jonas 

Technical  Problem 

The  goal  of  this  research  was  to  improve  our  understanding  of  the 
dynamic  structure  of  liquids,  disordered  solids  and  polvmers  by  studying 
motion  and  interactions  at  a molecular  level.  Simple  polyatomic  molecular 
liquids  and  water  were  of  main  interest.  In  this  research  summary  we  em- 
phasize in  particular  the  work  and  results  on  dense  liquids  composed  of 
molecules  approaching  spherical  shape  because  these  studies  enabled  us  to 
draw  a number  of  conclusions  of  general  validity  about  the  dynamic  structure 
of  these  liquids.  In  a sense  our  experiments  on  disordered  organic  solids 
and  polymers  must  be  still  regarded  as  exploratory  and  in  fact  represent  an 
initial  successful  stage  of  our  systematic  effort  in  this  specific  direction. 

General  Methodology 

In  order  to  obtain  detailed  information  about  molecular  motions  in 
liquids  and  provide  meaningful  experimental  tests  for  current  theories  of 
liquids  one  has  to  separate  the  effects  of  temperature  and  density,  i.e.,  one 
lias  to  use  both  temperature  and  pressure  as  experimental  variables.  It  is 
not  surprising  to  find  that  the  volume  changes  often  play  a dominant  role  in 
affecting  molecular  motions  in  dense  liquids  because  of  the  close  packing  of 
the  molecules.  The  capability  of  carrying  the  NMR  and  Raman  experiments  at 
high  pressure  over  a wide  range  of  temperatures  represents  an  important  and 
unique  feature  of  our  experiments. 

Nuclear  magnetic  resonance  relaxation  experiments  provide  detailed 
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information  about  reorientational  motions,  angular  momentum  correlation 
times  and  self-diffusion  coefficients.  Since  the  NMR  experiments  yield 
the  zero  Fourier  transform  of  the  correlation  functions  for  the  appro- 
priate motions,  we  started  using  laser  Raman  scattering  experiments  to 
obtain  detailed  time  dependence  of  the  correlation  functions.  The  question 
of  applicability  of  various  hydrodynamic  equations  on  a molecular  level  was 
of  interest  and  therefore  shear  viscosities  were  also  measured  as  a function 


of  temperature  and  pressure. 


Technical  Results 

The  results  of  our  studies  provided  direct  experimental  evidence  on 
the  validity  of  the  idea  that  the  dynamic  structure  of  a dense  liquid  is 
dominated  by  the  harsh  repulsive  intermolecular  forces  and  that  the  slowly 
varying  attractive  forces  can  be  neglected.  The  experimental  data  were 
analyzed  in  terms  of  the  rough  hard  sphere  model  of  liquids  based  on  the 
modified  Enskog  theory.  A rough  hard  sphere  model  of  a liquid  is  a model 
for  spherical  particles  for  which  interparticle  binary  collisions  occur  in- 
stantaneously and  are  capable  of  changing  angular  as  well  as  linear  momentum 
of  a particle.  In  addition,  the  rotational  and  translational  motions  are 
coupled.  Our  experiments  established  that  the  rough  hard  sphere  model  pro- 
vides good  description  of  the  dynamics  in  dense  liquids.  The  results  provided 
experimental  verification  of  the  molecular  dynamics  calculations  for  hard 
spheres  as  reported  bv  Alder  et^  al_.  From  the  experimental  relaxation,  dif- 
fusion and  viscositv  data  we  obtained  effective  hard  sphere  diameters  and 
their  temperature  dependence  for  a variety  of  molecular  liquids  (e.g.,  benzene. 


tetramethylsilane,  sulfur  hexafluoride,  etc.). 
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emphasizing  was  the  experimental  finding  that  the  time  dependence  of  the 
angular  momentum  correlation  function  can  be  approximated  as  a series  of 
uncorrelated  binary  collisions.  In  contrast  our  results  showed  the  im- 
portance of  many  body  correlation  effects  on  diffusion  and  shear  viscositv 
in  agreement  with  the  results  of  the  molecular  dynamics  calculations.  We 
showed  the  limitations  of  the  hard  sphere  model  at  high  packing  where  pre- 
dictions from  this  simple  model  strongly  deviated  from  experimental  data 
for  all  molecular  liquids  studied.  The  hydrodynamic  Stokes-Einstein  equa- 
tion was  found  valid  in  the  slipping  boundary  condition  for  the  liquids 
investigated . 

The  experiments  on  liquid  water  and  heavy  water  showed  that  compres- 
sion leads  to  significant  distortion  of  the  hydrogen  bond  network  with  the 
important  result  that  the  dynamic  behavior  of  water  under  high  compression 
resembles  that  of  normal  molecular  liquid  of  comparable  molecular  size.  At 
high  densities,  the  hard  core  repulsive  interactions  begin  to  dominate  over 
the  directional  interactions  which  are  mainly  responsible  for  the  open 
structure  of  liquid  water  at  low  temperatures  and  pressures. 

The  results  of  two  studies  were  published  in  the  area  of  disordered 
organic  solids  and  elastomers.  The  self-diffusion  study  of  the  disordered 
fee  crystalline  phase  of  neopentane  indicated  that  self-diffusion  proceeds 
via  a relaxed  vacancy  mechanism.  The  pressure  effects  on  the  glass  transition 
were  of  main  interest  in  the  NHR  study  of  several  elastomers  including  natural 
rubber. 

For  clarity's  sake,  the  detailed  discussion  of  technical  results  is 


divided  into  three  sections.  The  first  section  deals  with  the  experimental 
and  theoretical  results  obtained  in  our  studies  of  dense  liquids  composed  of 


39 


s' 


molecules  approaching  spherical  shape.  This  part  of  our  discussion  can  be 
regarded  as  a final  report  on  this  specific  class  of  liquids  because  we 
arrived  at  a number  of  conclusions  of  general  validity  and  further  progress 
in  this  area  must  come  from  development  of  better  theoretical  models.  Tin- 
second  section  is  devoted  to  hydrogen  bonded  water  and  heavy  water  which 
exhibit  many  anomalous  properties  when  compared  to  normal  atomic  and/or 
molecular  liquids.  The  third  section  summarizes  the  results  on  disordered 
solids  and  polymers  and  serves  mainly  as  an  illustration  of  future  directions 
of  our  studies 

i ) Molecular  Liquids 

The  NMR  spin-lattice  relaxation  times,  self-diffusion  coefficients, 
shear  viscosities  and  PVT  data  have  been  measured  in  a varietv  of  simple 
polyatomic  molecular  liquids  over  a wide  range  of  temperatures  and  pres- 
sures. The  dense  liquids^-^  such  as  benzene,  pyridine,  tetramethyls ilane, 
neopentane,  perf luorocyclobutane , sulfur  hexafluoride,  carbon  tetraf luoride, 
etc.,  are  composed  of  molecules  close  to  a spherical  shape.  The  analysis  of 
tlie  NMR  relaxation  data  and  transport  data  was  based  on  molecular  dynamics 
calculations  of  transport  coefficients  for  a hard  sphere  fluid  as  performed 
by  Alder  e_t  al^. 7 and  on  the  rough  hard  sphere  model  of  a liquid  as  proposed 
by  Chandler. R These  models  are  based  on  the  Enskog  theory g and  the  ideal 
that  in  dense  liquid  the  molecules  are  so  closely  packed  that  the  dynamic 
structure  is  dominated  by  the  intermolecular  repulsive  forces  depending  on 
the  shape  of  the  molecules.  As  far  as  our  results  are  concerned  the  pre- 
dictions of  the  rough  hard  sphere  model  agree  well  with  the  experimental  data. 
Our  data  confirm  the  molecular  dynamics  calculations  for  diffusion  and  shear 
viscosity  of  a hard  sphere  fluid  and  clearly  show  that  many  body  correlations 
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affect  in  a major  way  the  self-diffusion  coefficients  and  cause  deviations 
from  predictions  based  on  the  simple  Enskog  theory.  At  intermediate  den- 
sities the  self-diffusion  coefficient  of  a liquid  is  enhanced  due  to  long 
persistence  of  velocity  currents  (vortex  effect),  whereas  at  high  densities 
the  self-diffusion  coefficient  is  much  lower  than  the  Enskog  self-diffusion 
coefficient.  This  effect  can  be  ascribed  to  backscatteri ng  caused  by  the 
high  probability  of  reversal  of  velocity  of  a particle  upon  collision  with 
its  nearest  neighbor  in  dense  systems.  The  experimental  data  enabled  us  to 
calculate  the  effective  hard  sphere  diameters  for  the  various  liquids  studied 
and  also  obtain  their  temperature  dependence.  For  all  liquids  we  found  that 
at  high  packing  the  predictions  based  on  hard  sphere  model  deviate  signifi- 
cantly from  the  observed  values  of  self-diffusion  coefficient  and  shear  vis- 
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cosities. 

There  are  several  possible  reasons  for  the  interesting  behavior  of 
the  self-diffusion  coefficients  at  high  densities.  First,  there  is  a dis- 
tinct possibility  that  the  extrapolated  values  from  the  molecular  dynamics 
calculations  may  be  in  serious  error  at  high  densities.  Secondly,  the  ef- 
fective hard-core  diameters  can  become  a function  of  density  at  the  high 
packing  fractions.  This  possible  dependence  of  o upon  density  may  have  its 
origin  in  the  following.  One  has  to  realize  that  the  rough-hard  sphere  model 
is  still  a qualitative  model,  and  that  the  real  liquids  we  deal  with  in  this 
study  could  perhaps  be  represented  as  the  soft  sphere  liquid.  A softness  of 
the  potential  may  have  a significant  effect  on  the  self-diffusion  coefficients 
because  the  relationship  between  the  hard  spheres  and  the  soft  spheres  is  not 
well  understood  for  the  transport  processes.10 


In  contrast  to  the  behavior  of  the  self-diffusion  coefficient  and  shear 
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viscosity  no  correlation  effects  were  observed  for  the  angular  momentum 
correlation  t imes . 2 ’ c> ' 6 * 1 1 The  angular  momentum  correlation  time,  t , re- 
flects the  time  for  a molecule  to  lose  memory  of  its  initial  angular  mo- 
mentum and  as  such  it  may  be  expected  to  be  closely  related  to  the  molecular 
collision  frequency.  These  correlation  times  were  obtained  from  NMR  measure- 
ments for  liquid  systems  where  spin-rotat i on  interactions  dominate  the  relaxa- 
tion mechanism.  For  all  liquids  studied  (SF  , C.F  , CF, , CFCL.  C,DcF)  we 
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found  that  the  experimental  data  can  be  analyzed  in  terms  of  the  rough  hard 
sphere  theory  assuming  that  the  process  determining  the  time  dependence  of  the 
angular  momentum  correlation  function  can  be  approximated  as  a series  of  un- 
correlated binary  collisions.  Tt  is  satisfying  to  note  that  the  hard  sphere 
diameters  obtained  from  the  analysis  of  spin- rotat ion  relaxation  data  agree 
well  with  those  obtained  from  diffusion  and  viscosity  experiments.  One  should 
emphasize  that  in  our  studies  we  provided  a direct  experimental  evidence  that 
the  effective  hard  sphere  diameter  is  temperature  dependent.  This,  of  course, 
makes  physical  sense  because  the  repulsive  part  of  the  intermolecular  potential 
in  a real  fluid  has  a finite  slope  and  molecules  with  higher  kinetic  energy 
come  closer  to  each  other  during  a harsh  collision  and  thus  exhibit  a smaller 
hard  sphere  diameter. 

Since  our  experiments  provided  self-diffusion  coefficients  and  shear 
viscosities  of  various  liquids  over  a wide  range  of  temperatures  and  densities 
it  was  possible  to  test  the  validity  of  the  hydrodynamic  Stokos-F.instein  equa- 
tion on  a molecular  scale. This  relation  states  that  the  product  of  the 
diffusion  coefficient  and  the  viscosity  is  a constant  related  to  the  diameter 
of  the  spherical  particle  moving  through  a fluid.  The  many  body  correlations 
which  led  to  a decrease  of  diffusion  below  the  F.nskog  theory  are  compensated 


for  by  correlations  which  enhance  viscosity  by  nearly  that  much.  Over  a 
wide  range  in  volume  and  temperature  the  Stokes-Einstein  equation  is  re- 
markably precise  and  shows  that  the  slipping  boundary  condition  is  satis- 
fied well  for  most  of  the  molecular  liquids  studied.  It  is  interesting  to 
note  that  this  is  true  in  spite  of  the  fact  that  the  velocity  autocorrela- 
tion function  is  not  exponential  as  required  by  the  application  of  the 
Stokes  theory. 

i i ) Water  and  Heavy  Water 

It  is  well  known  that  water  exhibits  anomalous  behavior1^  with  pres- 
sure when  compared  with  normal  atomic  or  molecular  liquids.  All  the  anom- 
alies and  the  consequent  "strange"  behavior  of  water  under  compression  have 
their  origin  in  the  open  tetrahedral  structure  of  water  at  lower  temperatures 
due  to  a well  developed  hydrogen  bond  network.  We  measured1  _ 1 ’’  the  proton 
and  deuteron  spin-lattice  relaxation  times,  self-diffusion  coefficients, 
shear  viscosity  of  water  and  heavy  water  o>'er  a wide  range  of  pressures  and 
temperatures.  Thanks  to  the  capability  of  performing  our  experiments  under 
high  pressure  we  were  able  to  take  advantage  of  the  phase  diagram  of  water 
and  heavy  water  and  carry  out  experiments1'1  on  liquid  water  at  temperatures 
below  0°C  down  to  -15°C.  In  this  region  the  extremes  of  our  measurement  were 
ice  I and  a high  pressure  ice  V. 

Analysis  of  our  experimental  NMR  relaxation,  diffusion  and  viscosity 
data  lead  to  several  important  conclusions  about  the  motional  behavior  of 
water  molecules  in  compressed  liquid  water  and  heavy  water.  First,  the  initial 
increase  in  density  at  low  temperatures  produces  a faster  motion  of  the  water 
molecules.  This  is  easily  understood  in  terms  of  a simple  physical  picture. 
Since  reorientation  and  diffusion  of  water  molecules  must  proceed  via  breaking 
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and  reforming  of  hydrogen  bonds,  one  can  readily  see  that  high  compression 
which  changes  the  optimal  tetrahedral  hydrogen  bond  network  will  facilitate 
such  dynamic  processes.  Only  at  higher  densities  or  higher  temperatures 
further  compression  will  hinder  molecular  motions  due  to  increased  packing 
and  stronger  short-ranged  repulsive  interactions.  Under  these  conditions 
tlie  behavior  is  already  analogous  to  that  found  for  normal  atomic  and  mo- 
lecular liquids  where  compression  slows  down  the  motional  processes.  Second, 
compression  results  in  a decreased  coupling  between  the  rotational  and  trans- 
lational motions  as  reflected  in  the  relationship  between  the  reorientational 
correlation  time  of  water  molecules  and  shear  viscosity.  Therefore,  the  Debye 
equation  fails  to  describe  the  effect  of  density  on  the  reor ientat ion  of 
water  molecules.  Third,  we  find  that  the  activation  energies  for  relaxation, 
diffusion  and  shear  viscosity  decrease  with  increasing  density.  Fourth,  at 
high  temperatures  and  high  compression  the  dynamic  structure  of  water  re- 
sembles that  of  a hard  sphere  system.  From  these  findings  we  conclude  that 
compression  of  water  produces  ma-jor  changes  in  the  random  hydrogen  bond  net- 
work with  the  result  that  water  begins  to  behave  more  as  a normal  liquid 
under  high  compression.  All  our  experimental  results  are  in  agreement  with 
recent  molecular  dynamics  calculations  on  water  by  Still inger  and  Rahman.1' 

iii)  Disordered  Solids  and  Polymers 

Due  to  its  spherical  shape  neopentane  (C(CH^)^)  represents  an  excellent 
example  of  a molecule  exhibiting  a high  temperature  disordered  solid  phase,  a 
so-called  plastic  phase.  The  dominant  feature  of  this  fee  lattice  is  the 
large  degree  of  reorientational  freedom  but  there  is  also  evidence  suggesting 
a substantial  amount  of  translational  freedom.  The  NMR  proton  spin-lattice 
relaxation  times  have  been  measured  in  the  disordered  crystalline  phase  of 
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neopentane17  as  a function  of  temperature  and  pressure.  The  molecular  re- 
orientation and  self-diffusion  contributions  to  observed  relaxation  time 
were  separated  and  the  activation  parameters  were  calculated  for  these  two 
processes.  In  contrast  to  weak  pressure  dependence  of  the  overall  molecular 
reorientation,  the  self-diffusion  exhibited  a strong  pressure  dependence. 

The  activation  volume  for  self-diffusion  was  found  to  be  close  to  the  mo- 
lecular volume  of  neopentane.  Based  on  analogy  with  close-packed  metals 
this  finding  indicates  that  self-dif fusion  in  the  plastic  phase  of  neopentane 
proceeds  via  vacancy  mechanism. 

The  work  dealing  with  the  high  pressure  NMR  study  of  molecular  motions 
and  glass  transition  in  several  elastomers18  can  be  regarded  as  a successful, 
preliminary  study  for  future  experiments  related  to  glass  transition  in 
amorphous  polymers.  The  NMR  proton  spin-lattice  relaxation  times,  Tj , have 
been  measured  in  natural  rubber,  synthetic  cis-l,4-polyisoprene,  cis-1 ,4- 
pol vbutadiene  and  natural  rubber  gum  stock  over  a range  of  pressures  from 
1 bar  to  3 kbar  at  temperatures  from  -10°  to  100°C  at  24.2  and  56.4  MHz. 

Using  the  concept  of  the  most  probable  correlation  times  for  the  motion  of 
main  chain  segments  the  activation  enthalpy  and  activation  volume  were  cal- 
culated from  the  temperature  and  pressure  dependencies  of  the  T^  minimum. 

The  finding  that  AV4  increases  with  increasing  temperature,  and  the  relative 
role  of  thermal  and  volume  effects  on  relaxation  were  briefly  discussed.  The 
observed  pressure  dependence  of  the  glass  transition  temperature  for  these 
elastomers  indicated  that  configurational  entropy  rather  than  volume  plays  an 
important  role  in  determining  the  glass  transition. 

As  it  is  clear  from  the  above  discussion,  our  NMR  relaxation  experi- 
ments were  performed  under  conditions  of  extreme  pressure  and  temperature. 
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Changes  in  those  conditions  produce  subsequent  changes  of  electrical  pa- 
rameters in  the  NMR  sample  probe,  which  in  turn  lead  to  variations  of 
effective  radio-frequency  field  and  its  phase.  Thus,  frequent  reoptimiza- 
tion of  rf  pulse  lengths  and  phases  are  always  necessary.  Therefore  we 
decided  to  develop  a computer-directed  automated  system1  1 for  the  measure- 
ment of  the  NMR  snin-lattice  relaxation  times.  This  system  automatically 
generates  and  optimizes  the  pulse  sequences  and  carries  out  the  measurement 
of  spin-lattice  relaxation  times  automatically,  "e  found  three  major  ad- 
vantages of  this  system:  i)  substantially  improved  precision;  ii)  reduction 

of  measurement  time;  iii)  more  accurate  optimization  of  pulse  lengths  and 
phases  than  for  manual  operation.  Further  work  is  in  progress  on  an  auto- 
mated system  for  NMR  relaxation  experiments  in  the  rotating  coordinate  frame. 
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Final  Project  Report 
Structure  of  Liquids 
W.  H.  Flygare 


Technical  Problem 

j 

The  purpose  of  this  research  is  to  examine  the  dynamics  and  static- 
structure  of  liquids  including  liquid  crystal  systems  above  the  liquid  crystal 
phase  transition  temperature.  We  are  primarily  interested  in  the  inter- 
molecular  forces  which  are  important  in  determining  the  orientational  structure 
in  liquids.  However,  we  have  also  examined  the  behavior  of  a binary  system 
with  respect  to  hard  core  repulsive  forces  and  their  role  in  dynamic  prop- 
erties such  as  mutual  diffusion. 
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General  Methodology 

The  general  experimental  technique  is  light  scattering.  Starting  with 
a monocnromat ic  optical  source  (Ar+  ion  laser),  both  the  intensity  and  spectrum 
of  the  scattered  light  are  examined  to  gain  information  about  the  static,  and 
dynamic  features  of  scattering  medium.  The  above  data  is  interpreted  by  using 
a correlation  function  theory  for  the  scattering  medium.  The  experimental 
correlation  function,  which  is  the  Fourier  transform  of  the  scattered  light 
spectrum,  is  written  (according  to  the  ergodic  hypothesis)  as  a function  of 
the  space-time-orientation  correlation  function  of  the  particles  in  the  scat- 
tering medium.  A model  of  the  scattering  medium  then  allows  a comparison  with 
the  experimental  data.  By  carefully  studying  the  correlation  function  as  a 
function  of  temperature  and  density  (in  a high  pressure  bomb)  we  are  able  to 
infer  the  nature  of  the  close  range  forces  causing  both  t lie  static  and  dynamic 
features  of  the  liquid  structure. 
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Technical  Results 

Building  on  our  earlier  work  on  p-methoxybenzlidene-n-butylanil i ne 
(MBBA)!  which  concluded  that  orientational  pair  correlations  (OPC)  were  de- 
termined by  the  size  and  shape  of  the  molecule,  we  examined  several  chem- 
ical analogs  of  MBBA  which  changed  considerably  the  electric  dipole  moments 
in  these  molecules.*’  The  OPC  were  analyzed  as  a function  of  density  and 
temperature  as  a function  of  size,  shape,  and  moment  of  the  molecules.  We 
found  that  the  molecular  size  and  shape  (repulsive)  largely  determined  the 
OPC  in  these  liquid  crystal  forming  molecules. 

At  this  time  we  initiated  a study  of  intermolecul ar  forces  in  the  sys- 
tem of  substituted  benzenes  where  the  size  and  shape  remained  relatively  con- 
stant with  substitution,  but  the  electrical  moments  changed  markedly.  Before 
this,  however,  three  important  studies  were  necessary  in  interpreting  accu- 
rately the  OPC  data  in  these  systems.  These  studies  included  the  first 
quantitative  attempt  to  understand  on  an  experimental  and  theoretical  basis 
the  nature  of  the  "local  field"  at  the  scattering  molecule,1  a careful  study 
of  mutual  diffusion  in  binary  systems,  which  showed  for  the  first  time  that 
hard  core  diameters  could  be  transferred  from  one  binary  system  to  another, 
and  finally  a careful  study  of  the  polarizability  anisotropies  and  their  dis- 
persions was  also  completed.5  These  above  studies  were  a prelude  to  the  main 
study  of  orientational  pair  correlations  in  systems  of  equal  size  and  shape 
with  different  electric  dipole  moments. c’ 

The  nature  of  the  intermolecular  forces  which  are  important  in  deter- 
mining the  magnitude  of  orientational  pair  correlations  (OPC)  in  liquids  was 
studied  by  combining  high  pressure  and  depolarized  Rayleigh  light  scattering 
techniques.  For  molecules  of  similar  shape  (toluene,  bromobenzene , and  nitro- 


benzene) we  find  that  the  magnitude  of  the  OPC  as  measured  hy  depolarized 
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light  scattering  is  not  simply  a function  of  the  packing  of  the  molecules. 

For  sixteen  monosubstituted  derivatives  of  benzene  of  similar  size  and  shape, 
we  find  that  the  magnitude  of  the  OPC  correlates  well  with  the  square  of  the 
molecular  dipole  moment.  This  correlation  of  OPC  with  dipole  moment  is  under- 
stood by  relating  both  the  molecular  dipole  moment  and  the  magnitude  of  the 
OPC  to  the  detailed  distribution  of  molecular  electronic  charge.  The  light 
scattering  reorientational  relaxation  times  for  toluene,  bromobenzene , and 
nitrobenzene  as  a function  of  pressure  have  also  been  measured  and  found  to 
be  linear  functions  of  viscosity. 
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Implications  for  Further  Research 

The  most  obvious  implication  for  future  research  would  be  to  try  and 
determine  the  effects  of  an  orientating  electric  field  on  the  orientational 
pair  correlations  in  liquids  composed  of  molecules  with  large  dipole  moments 
These  studies  may  have  some  bearing  on  the  nature  of  molecular  pyroelectric 


materials  such  as  polyvinylidene  difluoride.  Can  molecular  polarization  be 
frozen  into  an  amorphous  material  by  fast  cooling  in  an  electric  field? 
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